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The Josephson current between d-wave superconductors is investigated in the framework of tight- 
binding lattice model. The junction is modelled by a small number of connecting bonds. It is 
obtained that the Josephson current through one bond vanishes when at least one of the super- 
conductors has (110) interface-to-crystal orientation. Interference between the nearest bonds is 
appeared to be very important. In particular, it is the interference term that leads to the nonzero 
Josephson current for (110) orientation. Also, interference of two connecting bonds manifests itself 
in non-monotonic behavior of the critical Josephson current in dependence on the distance between 
the bonds. 
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Electronic transport through the junctions between 
high-T c cuprate superconductors has been the object of 
interest for many years. In particular, the dc Josephson 
effect has been intensively studied theoretically«i*2iMiit. 
These investigations were inspired by the fact that su- 
perconducting order parameter dominantly exhibits d- 
wave symmetry in high-T c cuprates. This was estab- 
lished by SQUID-like experiments^ and the tricrystal 
experiment^. Theoretical investigations predicted a num- 
ber of results, which are the consequences of d-wave na- 
ture of superconducting order parameter. For example, 
the zero-energy Andreev bound state (ZES) is formed 
at junction interfaces^, what gives rise to anomalous en- 
hancement of the critical Josephson current in d-wave 
superconductor junctions at low temperaturesiSi^. Also, 
for a mirror junction, where the order parameter on the 
both sides of the junction is rotated by the same angle 
in opposite directions, a non-monotonic temperature de- 
pendence of the critical current was predicted!^' 3 . The 
temperature dependence of the critical current was stud- 
ied experimentally on the grain boundaries with well- 
defined lattice orientatior&iSiiiiiS. For mirror junctions 
the non-monotonic behavior was found irA whereas in 
other cases a monotonic behavior was reported&i&ii. 
In addition, for the 45° asymmetric junction a sin 2(p- 
like current-phase relation was predicted^. However, in 
the experiments on the asymmetric junctions not all 
samples exhibit the predicted sin2(y9-like current-phase 
relationiZiiiii 

The theoretical investigations of the Josephson cur- 
rent for the junctions between d-wave superconductors 
were carried out on the basic of continuous approachi^ 
as well as making use of tight-binding lattice model 5 . 
The main results of these methods are consistent with 
each other for the case of planar junctions. However, the 
lattice-model approach can not only give a more realistic 
description of the electronic structure of the copper oxide 
planes of high-T c superconductors and allows to mimic 
the corresponding Fermi surfaces, but it also gives the 
possibility to study electronic transport through quan- 
tum point contacts of various types. The recent advances 



in the fabrication of nanoscale devicesi£*i2*i£ii& (for a re- 
view see—) has provoked a renewed interest in the de- 
tailed analysis of models involving a few conducting chan- 
nels. Josephson current through quantum point contacts 
has been investigated in a number of papers since the pi- 
oneer work by Beenakker and van Houten- . The theory 
describing a single-mode quantum point contact between 
two s-wave superconductors in a site representation has 
been developed inSLSSiS. In particular, inS single-mode 
quantum point contact is modelled by the only bond con- 
necting two s-wave superconductors. At the same time, 
to the best of my knowledge, the junctions between d- 
wave superconductors through a few connecting bonds 
have not been considered yet. Also, interference of the 
connecting bonds has not been investigated by now. 

The present paper is devoted to these issues. I consider 
two d-wave superconductors in a mean-field site represen- 
tation connected by the only bond or by several bonds in 
the tunnel limit. It is shown that for the case of several 
connecting bonds their interference is very important and 
leads to the non-monotonic (oscillating) behavior of the 
critical current in dependence on the distance between 
the bonds. Furthermore, it is found that for the (110) 
interface-to-crystal orientation of at least one of the su- 
perconductors the Josephson current through each sepa- 
rate bond vanishes due to the symmetry and the current 
is entirely determined by the interference term. 

The principal scheme of the junction under consider- 
ation is shown on the Fig^ Two superconductors are 
connected by the N bonds with appropriate hopping ele- 
ments ii,...,tj\r. Then the full Hamiltonian of the system 
takes the form 



Tt = Hi + Tt r 



V. 



(1) 



where 7i; iT . correspond to the separate half-spaces and V 
contains coupling between them. For the each separate 
half-space we can use the usual mean-field lattice Hamil- 
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FIG. 1: The scheme of the considered junction. Left and right 
superconductors Si and S r are connected by the N bonds with 
appropriate hopping elements ti,...,t/v 



tonian in the tight-binding model 

Hl,r = -^A*4( r )c CT (r) - 51 *r,r'4(r) c <r( r ') 



•X){^4(r) C J(r') + fc.c.}. (2) 



Here /i is the chemical potential; t\ r ^ r are the hopping 
elements. For d-wave superconductors pairing is assumed 
to be nonzero only for the nearest neighbors A r r ± a = 
— A r r -(-b = A . Here r - are the site positions; a, b - are 
the basic lattice vectors (|a| = |b| = 1). The coupling V 
between the superconductors is modelled by 

V = - E {tr.^A^A^ + h.c.}. (3) 

a,i=l...N 

As we are interested in the Josephson current in such a 
system, the phase difference <p between the superconduc- 
tors should be taken into account. It is convenient to 
include it into the coupling term V: 



f = f. P #/2 

r i.li r i.r — 1 ' 
f — f . p-»0/2 



(4) 



Then we can assume the superconducting gaps A^ r and 
hopping elements ti to be real numbers in the whole sys- 
tem. 

Let G(r,r') = G(x,x' ,u> m ) be the Green's function 
of the whole system; Go(r,r') is the Green's function of 
the uncoupled superconductors; Gfc(r,r') = Gfc(r — r') is 
the Green's function in the bulk of the superconductors. 
Then we can represent G via Go and T-matrix: 

G(r,r') = G (r,r')+^Go(r,r^)f;^Go(r ji ^r'). (5) 

<X,B=l,T 

i,j=l...N 

Here all elements are 2x2 matrices in particle-hole 
space. Then it can be derived making use of Eq.© that 



T-matrix obeys the following equation: 



, ., ^ G (r i>a ,r k „)TZf=t?5 id 5 &tl3 . (6) 



T = (,r 
k = l...N 



Here I — r, r = I, and matrices tf are defined by 



t~e i0/2 
-he-**' 2 



(7) 



By the definition the Green's function &q(t, r') corre- 
sponds to the uncoupled superconducting half-spaces. 
Then we can write 



G o (r ; ,r r ) = 0, G o (r r ,r ; ) = 0, 



(8) 



for any 17 € Si and r r E S r . After doing this, we can 
resolve Eq.JSJ) in (Z, r)-space: 



rpr,r jr<r,l 

G l (l-G r & )-H (l-G l G%)-H* 
(l-G r G l )-H G r (l - G l G r )-H* 



(9) 



Here symbols with the check are 2N x 2N matrices in 
particle-hole and i,j = 1...N spaces. Namely T a -^ = 
T?f, i = il6 itj and Gg = tf G (r i>a ,r jja ). 

Let the y-axis be along the surface in the (a, 6)-crystal 
plane and the x-axis be the normal to surface. Then 
the functions Gq i ^ = Go(r,-, a , rj i<x ) can also be obtained 

from Gb by the T-matrix technique 2 ^: 



G (x, x', ky) = G b (x - x', k y ) 



(10) 



G b (x - x 0l k yi uj) G b (0,k y ) G b (x - x' ,k y ). 

Here xq is the position of isolating barrier between the 
superconductors (see24 for details) and Gb(x,k y ) is the 
Fourier transform with respect to k x of the bulk Green's 
function G(,(k) (h= 1): 



7r/d 



G b {n, k y ) 



2tt 



G b (k x ,k y )e lk * xd dk x . (11) 



k x — — 7r/d 



I only consider two possible interface-to-crystal orienta- 
tions of superconductors: (100) and (110). Then d = 1 
for (100) orientation and d = 1/V2 for (110) orientation. 
Instead of the usual square Brillouin zone k a = [— 7T, it], 
kb = [— 7T, 7r] I now use the surface-adapted Brillouin zone 
given by k x = [— Tt/d, it/d] and k y — [— ird, ird]. Then 
for the following calculations we should transform the 
Green's function Gq(x, x' ', k y ) into the full coordinate rep- 
resentation Go(r, r') = Gq(x, x' , y — y'). 
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The Josephson current through the barrier can be ex- 
pressed by (G - is the upper left part of G, uj m — 
(2m+l)7rT) 

J(^) = -2ieTj2 (Wi.rG(iV,r M )- 

Ulm 1=1. ..iV 

*ri, r ,r M G ! (r| j i,ri )r )) . (12) 

Taking into account Eq.JSJ) it can be obtained from 
Eq.© that 

G(ri ir , r^i) = ijTjjGa j,i 

i =1 - N „ „, (13) 
G(ri,i,r iir ) = Yl G ijTj'jGn jj. 

3 = 1. ..N 

Then for calculating the Josephson current we only 
need and TT\ elements of T-matrix. 

Let us consider now the case of one connecting bond 
(i.e. N = 1). Then Ea. (|12|l can be easily writen explic- 
itly: 

J W .2,tV?^^, (14) 

where Fq V and Fq T are the off-diagonal elements of G 
and Gq in particle-hole space: 

Ac _ ( Gq Fq \ , . 

Uo - ^ F Q Gq" J • 

The denominator in Ea. (|14|l takes the form 

Z(t, 4>) = 1 - P(G£G + G r G l - FZFle* - 
FZFle-**) + i\GlGl - F5F5)(G l Q G l - F^). (16) 

This denominator leads to high-order powers of trans- 
parency (more then first order in D ~ |i| 2 ) and is respon- 
sible for the deviation of the Josephson current from the 
sinusoidal behavior J(<f>) ~ sin(^). But for the particular 
problem it is more important to consider the numerator 
of Ea. l|14|) . It includes the anomalous Green's functions 
of the coinciding space arguments Fff = Fo(r Q ,r Q ) and 
Fq (with a = I, r), for uncoupled superconductors. Here 
r l and r r - are left and right ends of the bond connecting 
two superconductors. But it is easy to obtain, that for d- 
wave superconductor with (110) (i.e. 45°) smooth surface 
these Green's functions are zero: F(r,r) = F(r,r) = 0. 
This takes place due to the symmetry of sites positions 
and appropriate hopping elements with respect to the 
reflection (y — j/,-) — > — (y — yt) near the considered site 
[pi, yi) and simultaneous changing sign of order parame- 
ter under the reflection. These symmetry relations result 
in impossibility of flowing the Josephson current through 
one-bond contact if at least one of the superconductors 
is a d-wave superconductor with (110) orientation. 

The vanishing of one-bond Josephson current between 
(110) d-superconductors is a consequence of d-wave sym- 
metry of order parameter and has no analogue for the 



junctions between s-wave superconductors, where the 
current is nonzero for one-bond contact. 

Any reasons which do not change the above symmetry 
of the system (for example, surface pair breaking) cannot 
change this statement. But if the symmetry is destroyed, 
the nonzero one-bond Josephson current arises. If asym- 
metry is small, then current is small also. The possible 
reasons, giving the nonzero Josephson current through 
one-bond contact between (110) d-wave superconductors 
are: 1) not smooth surface of the superconductor (ends 
of facets or surface roughness) at the distance of the or- 
der of £o from the bond between the superconductors; 

2) the impurities in the bulk or at the surface of the su- 
perconductor, also placed not far then £ from the bond; 

3) nonzero is-component of the order parameter or mag- 
netic field. 

Let us consider the tunnel limit, i.e. only take into ac- 
count the first order of the barrier transparency D ~ t 2 . 
This means that values ti are sufficiently small and we 
can neglect Gq in comparison with I in Eq. © . It is worth 
to note that this approximation fails at sufficiently low 
temperatures due to the divergence of the Green's func- 
tion Go at uj m — > if the particular surface-to-crystal 
orientation of d-wave superconductor leads to the forma- 
tion of zero-energy surface bound states. 

In the tunnel limit 

fr,l = fr t f l ' r =l (17) 

Then the expression for the Josephson current takes the 
form: 

J(0) = 2 lG T^ £ Uij (Fq ijH j,^**- 

w m i,j = l...N 

Fo jt i F o ■ (18) 

It can be seen that Ea. ljT%)l contains two physically dif- 
ferent parts. First of them is the sum of the terms, which 
are proportional to t\. It represents the simple algebraic 
sum of currents through each bond separately (compare 
with Ea.H14[)). But the other terms are proportional to 
titj and give the interference part of the current. For the 
case of the junction with (llO)-orientation the interfer- 
ence part is the only non-vanishing term. 

Now I turn to the case N = 2. This is the simplest 
case for studying of the effects caused by the interference 
between the connecting bonds. The results are presented 
for junctions between d-wave superconductors with (100) 
and (110) orientations. In following the superconducting 
order parameter is assumed to be spatially constant. Al- 
though surface pair breaking is large for (110) surface 
orientation, this simplification does not change my re- 
sults qualitatively. 

In order to calculate the Josephson current Eq. (|18fl we 
need the Green's functions F(x, x' , y—y 1 ) and F(x, x' , y— 
y ). They are obtained in the model of nearest neighbors 
with the parameters Ao(T = 0) = O.li, /u = 0.5t. For this 
set of parameters T c w 0.173i and A max (T = 0) « 0.35i. 
I present the results for T w 0.6T C , where A (0.6T C ) w 
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FIG. 2: The critical Josephson current in two-bond contact (in units of eAo(T = 0)/ft) plotted as a function of the distance L = 
nd y between the bonds. T = 0.6T C , where T c is the superconducting critical temperature, (a) Junction between superconductors 
with (100) orientation (b) The same for (110) orientation. 



0.9Ao(T = 0). The hopping parameters for the tunnel- 
ing bonds are taken to be equal 1 1 = ti = O.lt and the 
bonds are placed at the distance L — nd y from each 
other. Here d y = d^ 1 is the period of lattice along the 
surface. For (100) orientation d y — 1 and for (110) ori- 
entation d y = \/2. Bonds i connect the last surface lay- 
ers of sites of both superconductors, therefore we should 
take x = x' = xf 11 "^ or x = x' = x s r wr ^ , where xf 1 ^ 
are the ^-coordinates of the positions of the last layers 
in left/right superconductor. The value y — y' equals 
for non-interference terms and y — y' = ±L for interfer- 
ence terms. Then the critical Josephson current can be 
written as: 

J c = 8 e T?^(|Fo| 2 + |F L | 2 ). (19) 

Here the term with Fq = F(y — y' = 0) corresponds to 
the current through each bond separately, while the term 
with Fl = F(y — y' = L) corresponds to the interference 
part of the current. The critical Josephson current de- 
scribed by Ea. (|19fl is presented on Fig[3 As the current 
for (110) orientation is entirely due to the interference 
term, it goes to zero with increasing of the distance be- 
tween the bonds. At the same time for (100) orientation 
the current tends to the value determined by the sum of 
two independent bonds. The interference part of the cur- 



rent oscillates and decays with increasing L due to the 
vanishing \Fl \ at L — > oo. It can be seen from Fig |2| that 
the characteristic distance between the bonds to consider 
them to be independent is £o ~ t/A.Q. It is worth to note 
that the interference part of the critical current is always 
positive for both orientations considered. 

In conclusion, in this paper I have studied the Joseph- 
son current through a contact with small number of con- 
necting bonds between two d-wave superconductors. It is 
obtained that the Josephson current cannot flow through 
one-bond contact connecting d-wave superconductors of 
(110) surface-to-crystal orientation with smooth surfaces. 
The expression for the Josephson current in the junc- 
tion with arbitrary number of bonds in tunnel limit is 
found. It is shown that the interference between the 
nearest bonds is very important. The Josephson cur- 
rent in two-bond junction is calculated. It is obtained 
that the interference of two connecting bonds manifests 
itself in non-monotonic behavior of the critical Joseph- 
son current in dependence on the distance between the 
bonds and leads to the nonzero Josephson current for 
(110) orientation. 
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